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We investigate the effect of progressive geographical infiltration on the topology of trajectory net- 
works. Trajectory networks, a type of knitted network, are obtained by establishing paths between 
geographically distributed nodes while following an associated vector field. These systems offer 
tools for modeling adaptive growth, development, and pathology of biological, transportation, or 
communication networks. For instance, the nodes could correspond to neurons or axonal branching 
points along the cortical surface and the vector field could correspond to the gradient of neurotrophic 
factors, or the nodes could represent towns while the vector fields would be given by economical 
and/or geographical gradients. The geographical infiltrations correspond to the addition of new 
local connections between nearby existing nodes. As such, these infiltrations could be related to 
several real- world processes such as contaminations, diseases, attacks, parasites, etc. Combined with 
a mechanism for elimination of nodes and connections, infiltration can model growth, development 
and adaptive plasticity in neuronal networks. The progressive geographical infiltration effect is ex- 
pressed in terms of the degree, clustering coefficient, size of the largest component and the lengths 
of the existing chains measured along the infiltrations. We show that the maximum infiltration 
distance plays a critical role in the intensity of the induced topological changes. For large enough 
values of this parameter, the chains intrinsic to the trajectory networks undergo a collapse which is 
unrelated to the percolation of the network also implied by the infiltrations. (Copyright Luciano da 
F. Costa, 2008) 

PACS numbers: 89.75.Hc, 89.75.Fb, 89.75.-k 



'No one remembers what need or command or desire 
drove Zenobia's founders to give their city this form, 
which has perhaps grown through successive super- 
impositions from the first, now undecipherable plan. ' (I. 
Calvino, Inivisible Cities) 



I. INTRODUCTION 

Graphs and complex networks can be classified into 
two major categories: geographical and non- geographical 
ones. Whereas in the latter type of networks, nodes do 
not have specific positions, in the former, each node has a 
well-defined spatial position, expressible by respective co- 
ordinates. Several real-world networks are geographical 
in nature, including power distribution (e.g. [l|), tourism 
(e.g. 0), transportation (e.g. [H), biological networks 
(e.g. bone structure [H, gene expression expression 
and developing neuronal networks 0, [1] ) • They all share 
the property that, to various extents, spatial proximity 
between nodes plays a role in shaping the connectivity 
structure. Often in these networks, spatially close nodes 



have a larger probability of being connected. Sometimes 
the role of spatial position is more intricate. For in- 
stance, in neuronal network development, axonal path 
finding is directed by the cooperation of multiple fac- 
tors. These include mechanical ones, such as the pres- 
ence of a fissure, the expression gradient of molecules 
as positive, permissive, or negative guidance factors and 
the cell adhesion molecules involved in fasciculation 
In addition we need to consider 0] neurothropic factors 
that regulate neuronal survival, differentiation, and sig- 
naling 0, [HI H3 > the gradients of neurotransmission [l2| , 
and the interactions amongst these factors [III, 03] • To 
take these into account, dynamical vector representations 
need to be associated with network nodes, vertices, or 
geographical locations. We wish to incorporate these re- 
quirements into geographical networks. 

A variety of geographical networks have been proposed 
in the literature (e.g. [H, [H, QjJ El ) ■ A new family of 
networks , namely the knitted networks, was proposed re- 
cently [1{J, [2fJ to include all networks defined and com- 
posed by paths, i.e. sequences of edges without repetition 
of nodes. 
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In this article, we expand the family of knitted net- 
works by incorporating structures generated by trajecto- 
ries defining paths following a given vector field. More 
specifically, a set of nodes is distributed within a given 
domain (a 2D space in this article, but the extension to 
higher dimensions is immediate); one node is chosen as 
origin, and the respective trajectory (line of force) is ob- 
tained while the nodes which are closer than a given max- 
imum distance to the current point of the trajectory are 
sequentially incorporated into the path. This procedure 
is repeated several times, yielding a network with connec- 
tions aligned to the vector field. In other words, the paths 
correspond to approximations of the solutions of the dy- 
namical system represented by the vector field. Figure Q] 
illustrates two trajectory networks obtained from the vec- 
tor fields <p(x,y) = (y,x) and <j>(x,y) = (y, —x) (b). 

Trajectory networks represent a natural putative 
model for several real-world structures and phenomen, 
for instance neural growth and development, includ- 
ing axonal navigation 0, the establishment of neuronal 
connections under the influence of neurotrophic fields 
(e.g. [§, [l(| [HI ) , neurotransmitter diffusion [lj] and their 
relation with adaptive plasticity [14| , the growth of trans- 
portation systems under geographical and economical in- 
fluences (e.g. 'every path leads to Rome'), the growth of 
trees and roots under influence of trophic factors [2l| . 
the development of channel-based systems such as bone 
structure and the vascular system, amongst many other 
important systems. 

The focus of attention in the current work is to inves- 
tigate how the topology of trajectory networks, a geo- 
graphical type of knitted network, is affected as a con- 
sequence of progressive geographical infiltration. By ge- 
ographical infiltration (hence infiltration for short), wc 
mean any process which interconnects pairs of nodes. 
Infiltration affects several real-world systems, e.g. the 
appearance of cracks along channels, the establishment 
of new local routes between towns and cities, contamina- 
tions between vessels of fibers, gallery building by para- 
sites, intentional attacks, internal spreading of diseases, 
to cite just a few cases. In the current work, the infil- 
tration process is simulated by selecting nodes at ran- 
dom and connecting this node to all other nodes which 
are closer than a maximum distance D p . Therefore, the 
adopted infiltration corresponds to the progressive incor- 
poration of tufts of local connectivity. 

Here, we investigat the effects of progressive infiltration 
on the topology of trajectory networks by quantifying 
the degree, clustering coefficient, size of the largest com- 
ponent, as well as the number and length of the chains 
present in the network. A recent study highlighted chains 
as an important category of network motifs [22j. Real- 
world networks often contain several chains, in ways spe- 
cific to their structure and function. Thus, these net- 
works are possibly the first theoretical model to natu- 
rally incorporate these motifs. These motifs are a conse- 
quence of the linking of spatially distributed nodes along 
the trajectories defined by the given vector fields. The 



common trait in real-world network structure that these 
models represent particularly well is the presence of in- 
dependent paths, with relatively few collaterals. There- 
fore, it becomes particularly important to characterize 
the structure of trajectory networks before and after in- 
filtration by considering the number and length of the 
existing chains. Interestingly, the effect of infiltrations 
can be either bad or good, depending on each specific 
system. For instance, the incorporation of additional lo- 
cal routes is in principle beneficial for transportation and 
communication systems. On the other hand, the addition 
of local connections in biological networks (e.g. bone or 
neuronal networks) may have catastrophic consequences. 
Observe that in the latter situation the main purpose of 
the chains/fibers is actually to provide mutual isolation. 
In both cases, the quantification of the effects of the infil- 
tration over the topology of the respective networks can 
provide valuable information to be interpreted from the 
perspective of each problem. 

This article starts by presenting the basic concepts — 
including the generation of trajectory networks and the 
geographical infiltrations — and follows by describing the 
experiments and discussing the respectively obtained re- 
sults. 



II. BASIC CONCEPTS 

A complex network is a graph exhibiting a particularly 
intricate structure. The connectivity of a undirected, 
unweighted network can be completely represented in 
terms of the respective adjacency matrix K, such that 
each interconnection between two nodes i and j implies 
K(i,j) = K(j,i) = 1, with K(i,j) = K(j,i) = being 
otherwise imposed. The immediate neighbors of a node i 
are those nodes which receive an edge from i. The degree 
of a node i is equal to the number of its immediate neigh- 
bors. Two nodes are said to be adjacent if they share an 
edge; two edges are adjacent if they share one node. A 
sequence of adjacent edges is a walk. A path is a walk 
which never repeats a node or edge. The length of a walk 
(or path) is equal to the respective number of involved 
edges. The clustering coefficient of nodes i is calculated 
by dividing the number of interconnections between its 
immediate neighbors and the maximum possible number 
of connections which could be established between those 
neighbors. 

A connected component of a network is a subgraph such 
that each of its nodes can be reached from any of its other 
nodes |38j . A chain is a subgraph of a network such as 
that each of its nodes has degree 1 or 2 and not additional 
nodes of degree 1 or 2 are connected to it [12] • The length 
of a chain is given by its number of edges. Two measure- 
ments which can be used to characterize the chains in a 
given network include the number of such chains and av- 
erage and standard deviation of their respective lengths. 
Chains are naturally related to paths along the network. 
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III. TRAJECTORY NETWORKS 

A family of networks, namely the knitted complex net- 
works, was introduced recently [H, H3] incorporating all 
networks organized around the concept of paths. Two 
main types of knitted networks were initially identified: 
path-transformed and path-regular. The former subcate- 
gory of knitted complex network is obtained by perform- 
ing the start-path transformation 19] on a given network 
(star and path connectivities can be understood as duals, 
e.g. through the line-graph transformation). Therefore, 
networks with power-law distribution of path lengths can 
be obtained by star-path transforming Barabasi- Albert 
networks [23| . The second type of knitted complex net- 
works, namely the path-regular networks, is particularly 
simple and involves starting with a set of N isolated 
nodes and performing several paths encompassing all 
nodes. Path-regular networks have been found to ex- 
hibit marked similar properties between different config- 
urations or nodes in the same configuration (e.g. [2(| H3|)- 
An even more regular version of the path-regular net- 
work, with all nodes exhibiting identical degrees, was 
later reported in [2^, [26[ . 

Geographical networks are characterized by the fact 
that each of their nodes has a well-defined spatial po- 
sition. Geographical networks represent an important 
category of complex networks because several real-world 
structures are inherently embedded into 2D or 3D spaces, 
and their connectivities are strongly affected by prox- 
imity and spatial adjacency. Given a set of spatially 
distributed nodes embedded in a continuous space to 
which a vector field is associated, it is possible to ob- 
tain geographical networks whose connections are a con- 



sequence not only of the proximity between nodes, but 
also of the orientations implied by the respectively as- 
sociated vector field. Several real-world can be thought 
as involving a geographical distribution of nodes and as- 
sociated vector fields. For instance, the neurons along 
the cortical surface can be represented as a set of geo- 
graphically distributed nodes, while their connections are 
established to a great extent as a consequence of neu- 
rotrophic fields (e.g. electrical or chemical gradients). 
Systems of streets, roads and highways can also be un- 
derstood as involving a set of spatially distributed nodes 
(the intersections between routes), with the interconnec- 
tions being established in terms of the spatial proximity 
between nodes as well as geographical and economical 
fields (e.g. the trend to connect to a big city, to avoid 
a geographical obstacle or to follow level-sets of height). 
Several other natural and human-made complex systems 
can be modeled by trajectory networks. Trajectory net- 
works are related to gradient networks (e.g. [271 l28l [29j| ) . 
field interactions [1, [g, H3| , as well as dynamical systems 
(e.g. [HI, HH). In the present work, we understand tra- 
jectory networks as a particular case of knitted networks. 

The trajectory networks considered in the present ar- 
ticle are obtained as follows. First, a two-dimensional 
workspace of size L x L is defined, and a vector field 
4>{x,y) is associated to it. For simplicity's sake we as- 
sume that — L/2 < x,y < L/2. All networks considered 
henceforth in this work are obtained for the vector field 
<t>{x, y) = (y, x). N points are distributed along this space 
with uniform probability. A total of N p trajectories are 
then performed while obtaining each network. A start- 
ing point is randomly selected, and the respective line of 
force (always parallel to the vector field) is calculated by 



4 



using the Euler leapfrog numerical method (e.g. [HI). At 
each current time, if a new node is found at a distance 
not exceeding D p , that node is connected to the previ- 
ous node, and so on. As it is clear from the example of 
trajectory network shown in Figure [T] the combination 
of proximity and orientation constraints while perform- 
ing the connections yield networks incorporating several 
chains, which closely follow the vector field orientation. 
Different degrees of interconnectivity between and along 
the chains can be obtained by varying the total number of 
points and the parameter D p . Observe that the number 
of chains is is reduced for larger values of D p /N. Once 
all trajectories are performed, the isolated points can be 
removed (as adopted henceforth) or not (allowing further 
connections) . 



IV. GEOGRAPHICAL INFILTRATIONS 

Given a geographical network, several types of pertur- 
bations of its structure can arise as a specific consequence 
of its geographical nature, in the sense that nodes which 
are spatially closer may interfere one another. For in- 
stance, in a neuronal system, unwanted connections may 
appear between nearby neurons as a consequence of dis- 
eases. In transportation systems, it is only too natural to 
incorporate new local connections to the network. Sev- 
eral other types of geographical interferences are possible, 
including those arising as a consequence of contamina- 
tions, attacks, infiltrations, amongst many other. 

In this work we incorporate progressive infiltrations to 
a given network geographical network by selecting one of 
its nodes and connecting to it all other nodes which are 
not further than a maximum distance Di. 



V. RESULTS AND DISCUSSION 

A set of 30 trajectory networks was obtained for the 
field (j)(x,y) — (y,x). A total of 1000 nodes was ini- 
tially distributed within a squre region of side L = 100 
centered at (0,0), and N p = 100 trajectories were nu- 
merically calculated. Starting from a randomly chosen 
node, each node at a maximum distance D p — 2 from 
the current growing extremity of each trajectory was suc- 
cessively connected. An example of obtained trajectory 
network is shown in Figure [1] Each of the 30 networks 
underwent progressive infiltrations assuming Di = 5 and 
A = 10. Figure [3 shows four stages (100, 200, 300 and 
400) along the successive infiltrations for Di = 5. Ex- 
amples of the results of infiltrations with Di = 10 are 
depicted in figure [3] 

In order to characterize the alterations in the topology 
of the trajectory networks as they underwent progres- 
sive infiltrations, a set of measurements (e.g. [34j |) was 
taken along the process. These measurements included 
the average and standard deviation of the node degree, 



clustering coefficient, size of the largest connected com- 
ponent, and chain lengths along successive infiltration 
stages. Only chains longer than 3 edges were considered 
in the respective measurements. These chains were iden- 
tified by starting from each of the network nodes with 
degree 1 or 2 and following along both sides (in case of 
degree 2) until the respective extremities of the chains 
(nodes with degree 1 or larger than 2) were found (each 
detected chain was removed from the network in order to 
accelerate the processing of the remaining nodes). The 
results obtained for Di — 5 and Di = 10 are shown in 
Figure and [31 respectively. Figure [B] and [7J show the 
above measurements for all the 30 considered networks. 

It is clear from Figures 2] to [7] that, as could be ex- 
pected, the degree and clustering coefficient both in- 
creased as a consequence of the addition of the infiltration 
tufts. Both such increases are sublinear, with a steeper 
decrease in the rate of clustering coefficient increase ob- 
served for Di = 10 (Fig.[7j. The relative sizes of the max- 
imum connected components suffer an abrupt transition 
before the 160 first infiltrations (most of the transitions 
take place before that value) for both settings of Di, but 
is more abrupt for Di = 10. This change is related to the 
percolation of the chains in the original network. Another 
relatively abrupt change is observed for the path lengths, 
most of which stabilizing themselves at a value near 6 for 
Di — 5 and 4 for Di = 10. The interval from the start of 
the infiltrations until the average length of the chains sta- 
bilizes (as observed above) is called the period of collapse 
of the chains. Very few networks remained with large 
average chain lengths larger after 200 infiltrations. This 
confirms the fact, evident from Figure [3J that the tuft 
infiltrations tend to quickly eliminate most of the long 
chains in the trajectory networks (the chain collapse). 
For larger values of Di, after the collapse of the chains, 
the vector field influence on the network connectivity can 
be hardly distinguished by visual inspection, such as in 
Figures [3^b-d) . It is important to keep in mind that the 
fact that small values of Di tend to imply little effect 
over the chain structure of the trajectory networks is ul- 
timately related to the number N of initial nodes and the 
maximal distance D p considered for chaining the nodes 
during the construction of the networks. 

The two involved critical phenomena, namely the per- 
colation of the networks and the collapse of the chains, 
were investigated further in order to search for possible 
relationship between their respective onsets. In order to 
do so, transition points along the successive infiltrations 
were identified automatically. These points, respectively 
T p and T c , correspond to the first occurrence of the value 
1 for the relative size of the largest connected component 
and the first occurrence of the average chain length which 
is smaller or equal than 5, respectively. Figure [5] shows 
the respectively obtained distribution of T p and T c ob- 
tained for the 30 realizations of networks with Di = 10. 
It is clear from this figure that the two critical phenom- 
ena taking place in the considered trajectory networks 
seem to be largely independent, in the sense that no cor- 
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relation has been observed between their critical values. 
Interestingly, as shown in Figure El the collapse of the 
chains can take place before the respective percolation. 



VI. CONCLUDING REMARKS 

Geographical networks represent an important cate- 
gory of complex networks because of their natural po- 
tential for modeling a large number of real-world and 
human-made complex structures and systems. At the 
same time, the category of complex networks built up by 
paths, namely the knitted networks, constitutes an im- 
portant superclass of complex structures because of their 



intrinsic association with the concept of paths (as oppo- 
sited to star connectivity) and random walk dynamics 
(e.g. [H, USUI])- I n this work, trajectory networks have 
been understood to belong to the supercategory of knit- 
ted networks as a consequence of the fact that these struc- 
tures are the result of path generation processes. Tra- 
jectory networks constitute a special case, in which the 
paths tend to follow an associated vector field. Our main 
interest in the present work, however, consisted in inves- 
tigating how the topology of trajectory networks changed 
as a consequence of geographical infiltrations. While sev- 
eral types of attacks and perturbations have been consid- 
ered and investigated in complex network research, rela- 
tively less attention has been devoted to perturbations in- 
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(c) (d) 
FIG. 3: The network in Fig. Q] after 100 (a), 200 (b), 300 (c) and 400 (d) infiltrations with A = 10. 



trinsically related to geographical constraints, especially 
the adjacency and proximity between nodes. Yet, several 
important real- world and human-made systems are prone 
to this type of perturbations, ranging from the onset of 
unwanted neuronal connections to the incorporation of 
new local routes to transportation systems. 

The main contributions reported in this article are 
listed and reviewed in the following: 

Trajectory networks as a special case of knitted com- 
plex networks: We have defined trajectory networks as 
a novel sub-class of knitted networks. This type of geo- 
graphical knitted network corresponds to an interesting 
case where the connectivity is the consequence of both 
the proximity between nodes and the orientation of the 
underlying vector field. 



New type of perturbation of network structure: We con- 
sidered, possibly for the first time, perturbations (or 'at- 
tacks') to geographical networks which depend on the 
proximity between the spatially distributed nodes. We 
focused attention on 'tuft' infiltrations, where a node i 
is randomly chosen and all other nodes which are closer 
than a maximum distance Dj are connected to node i. 
This type of topological change can be related to sev- 
eral real-world effects such as unwanted neuronal tangles 
as a consequence of diseases, establishment of local con- 
nections in transportation networks, contaminations, and 
attacks. 

Qualitative changes resulting from infiltrations: The 
progressive infiltration of a trajectory network was in- 
vestigated in a systematic manner, considering 30 real- 
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(a) Degree, Di= 5 




(b) Clustering coeff., Di= 5 
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(c) Relative size max. conn, comp., Di= 5 




(d) Path lengths, Di= 5 
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FIG. 4: Measurements of degree, clustering coefficient, size of the largest connected component and chain lengths in terms of 
the number of infiltrations (identified as 'time') with Di — 5 for a network obtained for the vector field <j>(x,y) — (y,x). 



izations of networks obtained for the same configuration 
with respect to the vector field <p(x,y) — (y,x). The 
changes in the networks topology was monitored by tak- 
ing several measurements including the degree, cluster- 
ing coefficient, size of the largest connected component, 
as well as the particularly relevant lengths of the exist- 
ing chains. The latter measurements are especially im- 
portant because the trajectory networks are inherently 
composed by chains. While the degree and clustering 
coefficients underwent relatively smooth increases, the 
size of the largest component and average chain lengths 
were subjected to relatively abrupt changes related to 
the percolation of the network (in the case of the largest 
connected component) and to the collapse of the chain 
structure (in the case of the average chain lengths). The 
value of Di was found to be have great influence on such 
topological changes induced by the infiltrations, with val- 
ues much larger than D p implying particularly intense 
changes, especially regarding the chain structure. Af- 
ter the collapse of the chains, the effect of the original 
vector field on the network connectivity could hardly be 
discerned. Such findings are particularly important for a 



large number of real-world structures underlain by tra- 
jectory networks and geographical infiltrations. 

Independence of percolation and collapse of chains: 
The progressive infiltration of trajectory networks in- 
volves two critical phenomena: its percolation and the 
collapse of its chain structure. Interestingly, no clear re- 
lationship between these phenomena has been identified 
by considering the critical times T p and T c . This implies 
that the collapse of the chains can not be predicted from 
the percolation of the respective network, and vice- versa. 
As a matter of fact, it has also been observed that the col- 
lapse of the chains can take place before the percolation 
of the respective network. 

Modeling of Brain Development The current study il- 
lustrates that, in order to avoid pathological network 
conditions, besides growth, a mechanism for the selec- 
tive elimination of connections is also necessary. Such 
a mechanism can be observed at work in brain develop- 
ment. The formation of neuronal networks involves the 
extensive growth, but also elimination of neurons and 
connections (35|. Isolated nerve cells undergo apoptosis; 
dendritic arbors are being built and retracted based on 
signaling efficacy and electrical activity in the pre and 
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(a) Degree, Di= 10 




(b) Clus tering c o ef f . , Di= 1 
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(d) Path lengths, Di= 10 
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FIG. 5: Measurements of degree, clustering coefficient, size of the largest connected component and chain lengths in terms of 
the number of infiltrations (identified as 'time') with Di — 10 for a network obtained for the vector field <f>(x,y) = (y,x). 



postsynaptic neurons. As a result, synapses undergo ex- 
tensive rewiring after their initial attachment [3(| . These 
processes work together to maintain a functional network 
architecture for effective communication between brain 
cells [13. 

The several possibilities of future work include but are 
not limited to the following: 

Other types of vector fields: It would be interesting 
to investigate how the patterns of topological changes 
observed in this work extends to trajectory networks ob- 
tained by considering other vector fields, as well as other 
configurations of the involved parameters. 

Orthogonal infiltrations: In this work we focused at- 
tention on tuft infiltrations. It would be interesting to 
study the topological changes of trajectory networks with 
respect of other types of geographical perturbations, such 
as connecting points according to proximity and orienta- 
tions orthogonal to the vector field (possibly also through 
trajectories). 

Infiltration by increasing distances: While the infiltra- 
tions implemented in this article consisted in selecting 
nodes followed by tuft interconnection, it would be par- 
ticularly interesting to investigate the topological alter- 



ations of trajectory networks while all pairs of nodes are 
joined according to successive distances. Such a type 
of infiltration is guaranteed to completely eliminate the 
chains after a critical interval. 

Application to real-world networks: It would be inter- 
esting to quantify the alterations of real-world networks 
expressible by trajectory networks, including trans- 
portation networks, power distribution, communications, 
tourism and neuronal systems. 

Application to Image and Shape Analysis: The analy- 
sis of images containin g object s and shapes has remained 
a great challenge (e.g. [30, [37|). It would be particularly 
interesting to consider the application of the concepts 
and methods reported in the current work to such prob- 
lems. More specifically, trajectories can be obtained in 
gray-level images by considering their respective gradi- 
ent fields. So, by distributing points through the image 
and interconnecting them while taking in to account tra- 
jectories driven by the gradient fields, it is possible to 
obtain respective network representations incorporating 
a great deal of the intrinsic geometric features. Shapes 
represented by their contour can also be mapped into 
trajectory networks by considering vector fields induced 
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(a) Degree, Di= 5 




(b) Qnstering coeff., I>i= 5 
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FIG. 6: Averages of degree, clustering coefficient, size of the largest connected component and chain lengths in terms of the 
number of infiltrations (identified as 'time') with D; = 5 for each of the 30 networks obtained for the vector field <f>(x, y) = (y, x). 
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FIG. 8: Scatterplot of the percolation (T p ) and collapse (T c ) 
critical times for Di — 10. 



